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Abstract 

The model of the interaction of the noise from a spreading subsonic jet with a panel- 
stringer assembly is studied numerically in two dimensions. The radiation resulting from this 
flow/acoustic/structure coupling is computed and analyzed in both the time and frequency do- 
mains. The jet is initially excited by a pulse-like source inserted into the flow field. The pulse 
triggers instabilities associated with the invisc.id instability of the jet mean flow shear layer. 
These instabilities in turn generate sound which provides the primary loading for the panels. 
The resulting structural vibration and radiation depends strongly on their placement relative to 
the jet/nozzle configuration. Results are obtained for the panel responses as well as the trans- 
mitted and incident pressure. The effect, of the panels is to act as a narrow filter, converting the 
relatively broad band forcing, heavily influenced by jet instabilities, into radiation concentrated 
in narrow spectral bands. 
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1. Introduction 


In this paper we study numerically the interaction between jet noise, emanating from a high subsonic 
jet, with a flexible aircraft-type structure. The emphasis is on (1) jet noise generated by spatially 
growing instabilities associated with the shear layer of the jet as an excitation mechanism for the 
flexible structure, (2) the resulting vibration of and radiation from the flexible structure due to the 
flow/acoustic/structure coupling and (3) the effect of the placement of the panels, relative to the 
jet/nozzle configuration, on panel response and radiation. The precise details of the jet flow field 
are deemphasized via the use of a model jet flow field in the Euler equations. 

The numerical computations presented here are a step in determining the role of jet noise 
excitation on both structural fatigue (sonic fatigue) of aircraft panel assemblies and interior aircraft 
noise levels. The results demonstrate the important role of instability wave generated sound on 
the frequency response of and radiation from the ensuing structural vibration. In particular, since 
instability wave generated sound is highly directional, the panel placement relative to the jet/nozzle 
configuration is an important factor in determining the resulting panel response. 

The jet is excited by a pulse starter, an artificially imposed source that triggers continuous 
perturbations in the flow field. We concentrate on the long time behavior of the excited jet, i.e. the 
behavior of unsteady disturbances in the jet long after the initial disturbance generated by the pulse 
has propagated from the region of interest. The pulse starter is applied to the Euler equations, 
suitably modified so that only disturbances in a realistic spreading jet profile are computed. As a 
result of the pulse, an initial acoustic disturbance is generated. This acoustic disturbance propagates 
into the farfield while undergoing convection and dispersion due to the flow field. More importantly, 
the initial disturbance triggers the development of an instability wave, i.e., a large scale structure 
which is slowly convected downstream along the jet axis, due to interaction with the unstable jet 
shear layer. This instability wave, which is qualitatively similar to structures observed in real jets, 
is a source of jet noise as it triggers the development of other instability waves in a continuous 
fashion long after the initial disturbance due to the pulse has propagated away from the flow field 
and the panels. In particular, the qualitative features of large scale instability waves as sources of 
jet noise are simulated [7, 16, 22, 34]. The model also allows simulation of dispersive effects such 
as the frequency dependent refraction of sound by the flow field. 

The use of the pulse starter allows the qualitative simulation of instability wave generated 
sound in a spreading jet, a key component of real jet noise, with sufficiently reduced computational 
requirements so that high resolution computations of the fully coupled flow/acoustic/structure 
interaction can be readily obtained. This methodology is similar to that employed in [2, 21, 22] 
where the behavior of instability wave generated sound in subsonic axisymmetric jets (without any 
nearby structure) was considered. These results demonstrated that the starter pulse was effective as 
a computational tool to excite instability wave generated sound and to study the bending (refraction 
effect) of sound waves through the jet flow field. The results also demonstrated that the instability 
wave generated sound was highly directional, peaking at low to mid angles from the jet axis, and 
concentrated at low frequencies. 

There have been many important studies on jet noise. The sources of jet noise were recognized in 
the Lighthill acoustic analogy [17]. The resulting model led to scaling arguments for the dependence 
of the farfield sound on jet velocity, although the actual sources had to be modeled for calc ulation 
purposes. An independent theory of sound radiation, using a different formulation of the sources, 
was obtained in [29]. Extensions of this work to account for shear interaction terms and vortical 
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flows followed [18, 28, 33]. Experiments demonstrated the existence of large scale structures or 
instability waves in jets and the role of instability wave generated sound as an important component 
of jet noise [7, 16, 22, 34]. Analytical studies provided further elucidation of the role of instabilities 
as an important source of jet noise [5, 13, 14, 26, 27, 30] as did the computations [2, 21, 22]. 

More recent computations have included the effect of viscosity and turbulence models, for 
example, [25]. These computations demonstrated the substantial difficulties in simulating sound 
generation from spreading jets directly from the Navier-Stokes equations. In the approach taken 
here, the jet profile is modeled and included as a source term in the Euler equations. The resulting 
model allows the formation of inviscid instabilities and associated jet noise, which qualitatively 
models many of the observed features of the natural sources of sound in subsonic jets. The model 
also allows for the effect of convection and dispersion of acoustic waves due to the flow field. 

An understanding of the nature of the jet forcing is an essential prerequisite to an understanding 
of the role of jet noise on the structural vibration of aircraft panels. Since instability waves are 
one of the most important sources of jet noise, particularly for subsonic jets, it is important to 
determine the role of instability wave generated sound as a forcing mechanism for flexible aircraft 
panels. 

For sources that can be well approximated by modes of the panel, simple normal mode analysis 
suggests that the displacement would exhibit spectral peaks near the resonant frequencies of the 
panel, i.e., the natural frequencies. Even for such sources, other quantities such as the velocity 
and the radiated pressure, which must be obtained from solution of the wave equation (either fully 
coupled to the panel motion as done here or a posteriori using the uncoupled panel velocity as 
a boundary condition) will not necessarily exhibit spectral peaks at the natural frequencies. The 
spectrum of these quantities will depend on both the details of the forcing function and on the 
parameters of the panels. When the incident pressure is jet noise, as considered here, the situation 
is yet more complicated as the sources will be convecting downstream with the flow and may not 
be well approximated by modes of the panel. In this case even qualitative details of panel response 
would be expected to require specification of the sources accounting for the most physically relevant 
features of jet noise generation. 

Various approaches have been employed to determine panel response in a system exhibiting 
flow/acoustic/structure interaction. The sources can be taken from experimental measurements, 
(e.g., [8, 24]). Another approach is to model the sources using appropriate simplifications (eg., 
[6, 35]). In the former case, the panel response could then itself be measured experimentally as 
in [20] or computed via solution of the resulting panel equation from the measured sources as 
in [8, 24]. While recent experiments [20] in supersonic jets have succeeded in isolating features 
of panel response due to oscillating shocks and jet instability, in general the use of experimental 
measurements dictates that all physical mechanisms associated with the jet noise be intertwined, 
thus making it difficult to ascertain the role of individual physical mechanisms of jet noise on panel 
response. 

The modeling of sources necessarily involves significant simplifications of an extremely complex 
phenomena, namely the generation and propagation of jet noise, and also typically does not account 
for the full fluid/acoustic./structure coupling of the panel response. In the present work, we focus 
primarily on the role of sound generated from instability waves, large scale structures associated 
with the inviscid instability of the jet shear layer, as an excitation mechanism for the structural 
vibration. Such waves are known to be important sources of jet noise for both subsonic and 
supersonic jets [16]. The computations, while employing a model jet, are able to simulate the 
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qualitative features of such noise and the ensuing pane] response and radiation. We employ a 
computational model in which the panel response and resulting radiation are fully coupled to the 
evolution of disturbances in the jet flow field. Similar fully coupled calculations were employed 
in [9] where the response of a panel to a normally incident plane acoustic wave was considered. 
The computations simulated many of the features observed experimentally in [23]. The results 
demonstrated that highly complex responses, including subharmonic formation, quasbperiodic and 
possibly chaotic responses could be obtained for sufficiently high excitation level. 

Our results demonstrate that the forcing due to instability wave generated sound can lead to 
enhanced low frequency vibration and radiation from the panels. The effectiveness of jet instabilities 
in exciting the panel is highly dependent on the placement of the panel relative to the jet exit. The 
panel acts as a narrow filter in converting a relatively broadband forcing field, heavily influenced by 
jet instabilities, into a relatively discrete vibration and radiation spectrum. The resulting radiation 
is beamed primarily in the vertical and downstream directions. While the precise frequencies of 
the radiation will depend on the details of the flow field and the mechanism of exciting the jet 
shear layer instabilities, we believe that the general features of the panel response presented here 
illustrate the role of jet noise in exciting vibration and radiation from nearby surfaces. 

In section 2 we describe the model and the resulting numerical method. In section 3 we describe 
our results. In section 4 we summarize our conclusions. 

2. Problem Formulation 

The computational domain is shown in Figure 1. Fluctuating pressure, density and velocities are 
computed in two domains, that which contains the jet, exiting from a nozzle of width D, and the 
domain on the other side of the wall boundary. We will refer to the two domains as the jet and 
radiation domains respectively. The w'all boundary is a rigid wall containing two adjacent flexible 
cutout sections (denoted A and B in Figure 1) with rigidly clamped boundaries. The panels vibrate 
in response to excitation from jet noise and radiate sound into both domains. We focus primarily on 
acoustic radiation into the radiation domain, as the radiation into the jet domain is small compared 
to the large disturbances already p resent in the jet. 

The numerical method involves coupling the computation of a nonlinear equation governing the 
panel responses (the beam equation) to an Euler computation performed separately in both the 
jet and radiation domains. The difference in pressure across the panels (computed from the Euler 
equations) is the source term for the beam equation, and its solution provides the normal velocity 
to be used as a boundary equation for the Euler equations. We assume that the panel displacement 
is small relative to the acoustic wavelengths so that in each of the two Euler computations the 
normal velocity (obtained from the beam equation ) is imposed at y — 0. 

In the jet domain the Euler equations are modified to account for the jet flow. We assume 
a straight pipe of width D from which the jet exits. The solution is computed both within and 
exterior to the pipe, while rigid boundary conditions are assumed on each side. In the radiation 
domain, the Euler equations are directly computed. We note that for linear radiation, an alternative 
is to impose the radiation condition as in [32]. Thus, in the computations presented here, the jet 
flow field is fully coupled to the panel vibration and ensuing radiation. 

The nonlinear beam equation is given by 

c) 4 z () 2 Z d 2 Z dz + 

Dp M~ Nx W + pph di 2 + 1 Tt- p ~ p ' ( ) 
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where z represents the plate transverse deflection, p v the mass per unit volume of the plate, h the 
plate thickness, 7 the physical damping, and D p = Mh 3 / 12(1 — // 2 ) is the stiffness of the plate 
where M is the modulus of elasticity and v is the Poisson ratio of the plate material. The coefficient 
N x of the nonlinear term is given by 


N x 


2L J XQ \oxJ 


( 2 ) 


which represents the tension created by the stretching of the plate due to bending. The solution of 
(1) is obtained using an implicit finite difference method developed in [15] coupled with Newton’s 
method to account for the nonlinear term N x . Each panel is assumed clamped at both ends. 
The quantities p + and p ~ , the pressures above and below the panels, are determined from the 
computation of the radiation and jet domains respectively. 

The fluid equations are solved in conservation form. In terms of the vector 


w = (p,pu,pv, E) t 

the equations can be written in the form 


Wt + F x + G y — 0 (3) 

where subscripts denote derivatives, and the flux functions F and G are given by 

F = (pu,pu 2 + p,puv, u(E + p)) T , G = (pv,pvu,pv 2 + p, v(E + p)) T 

respectively. Here, p is the density, u, v are the x and y components of the velocity respectively 
and E is the total energy per unit volume, 

E = l/>(« 2 + v 2 ) + c-vpT , 

where T is the temperature and c v is the specific heat per unit volume. The pressure, p, is obtained 
from the equation of state 

p — pRT 

where R is the gas constant. We note that the equations are solved separately in each domain. We 
do not distinguish the vector w in each domain for clarity. 

We next restrict attention to the jet domain. Within this domain, the system (3) is modified 
to account for two different source terms 

Wt + F x + Gy = s-i(t,x,y) + s 2 (x,y). (4) 


The source term, 


h (t,x,y) = f f{t)g(x,y)w /, 


( 5 ) 


serves as the starter pulse to excite the jet. Here, e describes the amplitude of the source, f(t) is 
a pulse which describes the time dependence of the source, g{x y y) is a Gaussian peaked at a given 
location along the jet axis, and 


W/ = (1, Uo, 0, -p 00 ) T 

7 - 1 
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where the subscript “oo” denotes ambient quantities and Uo(x,y) is a model describing aspreading 
jet profile taken from [19] in Cartesian coordinates. The vector wj is chosen as follows. The compo- 
nent in the continuity equation serves as a mass injection. This is balanced by the component in the 
energy equation so that in the absence of flow and boundaries the solution remains homentropic. 
The component of the source in the axial momentum equation ensures that momentum is added 
to the system as the mass injected is convected downstream of the flow. The functional form of 
the pulse f(t) is given by 

f{t) = exp(-at' 2 - b/t 2 ), (6) 

with constants a and b. In the absence of flow and bounding surfaces, the solution would be a 
circular wave. As part of the validation of the numerical scheme, the computation with no flow has 
been extensively checked by grid refinement to give a convergent circular wave. 

The second source term y) is designed so that in the absence of the starter pulse, i.e. ( = 0, 
the solution to (4) would be a stationary profile corresponding to a spreading jet. The specified 
solution is 

Wo “ ( Poo ? Poo Uo i 0? Poo U 0 / 2 T C v T iyo ) 

In the limit of small the effect of the source term -s 2 leads to the Euler equations linearized 
about a spreading jet as in [21, 22] which by itself need not be a solution to the Euler equations. 
Incorporation of this term within the nonlinear Euler equations allows for the computation of 
nonlinear disturbances within the mean jet profile as in [2]. The inclusion of this source term 
separates the computation of the disturbance, in particular the resulting instability waves, from 
the computation of the mean flow (i.e., the spreading jet). Thus, the resulting system of equations 
allows for the simulation of instability waves and the resulting sound generation, together with the 
bending of acoustic waves in the jet flow field without requiring the computation of the spreading 
jet itself. Although this is a simplified model, the resulting system simulates many of the observed 
features of instability wave generated jet noise and permits high resolution computation of the 
coupling of jet noise with the flexible panels and the resulting radiation from the panels. 

The initial conditions are taken to be the mean state w 0 in the jet domain and ambient data 
in the radiation domain. The panels vibrate in response to excitation by sound from the jet as 
they would in an aircraft on the ground since the mean state corresponds to a static jet. Both 
the jet domain and the radiation domain are in principle unbounded regions. These domains are 
closed by the imposition of non-reflecting boundary conditions at the farfield boundaries, indicated 
in Figure 1. At these boundaries, we impose conditions designed to absorb outgoing circular waves 
[4]. Specifically we assume that the pressure has the functional form 

p(t,x,y) ~ hii-ovt. — r,9)jy/r , (7) 

where x and y are expressed in polar coordinates from a specified origin and c (X1 is the ambient 
sound speed. Differentiation of (7) with respect to t and r results in the relation 

Vt + '’ooPr + (p - Poc.)/(2r) = 0(r~ 3/2 ). ( s ) 

Neglecting the right hand side of (X) results in a radiation boundary condition which is effective 
in simulating outgoing waves provided the boundaries are sufficiently distant from all sources. We 
have verified that there are negligible boundary reflections for the data presented in this paper. 

The Euler equations are solved using an explicit finite difference scheme which is fourth order 
accurate in the spatial coordinates and second order accurate in time. The scheme is a generalization 



of the second order MacCormack scheme to allow higher order accuracy in space [10]. The numerical 
method has been validated on other problems of spatial instability. In particular, in [3] it was shown 
that the use of fourth order differencing can lead to a dramatic improvement in the calculation of 
growth rates for spatially unstable disturbances. In the numerical method we employ operator 
splitting in the two coordinate directions. The scheme is discussed in detail in other publications 
[ 21 , 22 ]. 

3. Results 

We consider a configuration as indicated in Figure 1. The jet exits from a straight nozzle of width 
D = 2in, assumed rigid on both the interior and exterior sides. The jet starter pulse (6) is located 
in the potential core at approximately 1.15 D downstream of the nozzle exit on the axis of the 
jet. The constants a and 6 in (6) are a = 3.4(c 00 /£)) 2 and b = 380 (c 00 /D)~ i , chosen to give a 
frequency spectrum for /' within the range of interest for subsonic jets. An infinite wall is located 
approximately 6 D above the jet and parallel to the nozzle. The wall is assumed rigid, except for 
two regions where flexible, aluminum, aircraft-type panels (panel A and panel B), with clamped 
boundaries are located. The panels are of length 5 D and thickness 0.01 D and are centered at 
x = 0 D and x = 5.24 D respectively. Other parameters of the panels are typical of aluminum. The 
origin of coordinates is chosen to be the horizontal location of the nozzle exit for x and the vertical 
location of the rigid wall for y. For the computations described below, we assume an exit flow Mach 
number of 0.65, with exit velocity Uj = 0.65 c^. Ambient conditions are assumed for p and p. Both 
the jet and radiation domain extend 48 D downstream from x = 0, 12 D in the upstream direction 
and 48Z) in the y direction. There are no detectable reflections from the artificial boundaries for 
the data presented here. We employ a grid of 351 x 351 points in the jet domain and 301 x 301 
points in the radiation domain. The grid in the jet domain is stretched to improve resolution of 
the jet shear layer and source region. The solutions described here have been validated by grid 
refinement. 

Our results are presented in three parts; the jet domain, including the flow and acoustic radiation 
from the jet, the responses of the panels and the acoustic radiation from the panels. 

3a. Jet Flow Domain 

Nonstationary behavior in the jet is triggered by the pulse starter, which generates a disturbance 
that propagates through the jet, interacts with the shear layer and then propagates into the farfield 
as sound. In the farfield this initial disturbance is non-circular due to convection. As a result 
of the interaction of the acoustic disturbance with the the shear layer gradient of the mean jet 
profile, instability waves are created. This was shown in [1], in the context of a linearization of 
(3). These waves are inviscid, a consequence of the fact that the Euler equations are employed in 
the computation. The instability waves propagate longitudinally along the jet at speeds less than 
the speed of sound, spreading transversally as does the mean flow. After an initial growth, their 
amplitude decays with the velocity as the jet spreads. As they propagate downstream they act as 
sources of additional disturbances as described in [17, 18, 26, 27, 29, 30] which propagate into the 
farfield as sound and also serve to excite further instability waves in the jet. This behavior is also 
consistent with experimental observations and with previous linear and nonlinear computations [2] 
and [22]. The resulting sound radiation forces the panels into a broadband response and radiation. 
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The acoustic radiation persists after the initial disturbance generated by the pulse starter has 
propagated a significant distance from the panels and away from the region of interest. 

Figure 2 shows Power Spectral Densities (PSDs) of the farfield pressure (36D from the nozzle 
exit) for the case of a domain without any wall at 30°, 60°, and 90°. In order to help clarify 
the effects of the flow, the figure also includes the results for a computation without any flow 
(but with a nozzle) at 90°. We note that the no flow computation is not quite omni-directional 
due to scattering from the nozzle. We have verified that in the absence of the nozzle the no flow 
computation is omni-directional to a high degree of accuracy. This figure illustrates the role of the 
jet as an amplifier of sound. The acoustic farfield is least affected by the flow for angles near 90°, 
a result also demonstrated in [22]. The data at 30° shows two peaks, one near 700Hz. and another 
near 220Hz. The 700 Hz. peak corresponds to the peak frequency of the pulse, amplified due to the 
presense of flow as indicated by the large difference with respect to the no flow case. The 220Hz. 
peak results primarily from the shear noise generated from flow instabilities. Such a behavior 
was also observed in [2, 22] for axisymmetric jets. Low frequency peaks are also observed at 60°, 
although the amplitudes are reduced with respect to the pronounced low frequency peak at 30°. 
This figure illustrates that the computed sound radiation has features consistent with experimental 
observation. 

The behavior of the resulting pressure field in the presence of the bounding wall, both within 
and outside of the flow, as well as the resulting acoustic radiation from the panels, is shown in 
Figures 3a-3e for successive values of time. Note that although both the jet and radiation domains 
in the computation extend 48 D in y y the figures show solutions extended only to 30 D in y in order 
to expand the visualization of the acoustic field. In interpreting these figures, note also that the 
actual levels assigned to each color for the jet and radiation domain differ so that different colors 
are assigned to the ambient regions. However, for each domain the values associated with the colors 
remain the same. 

Figure 3a shows pressure contours at t = 27.9 D/c^, indicating the initial development of the 
unsteady jet flow' and acoustic field of the jet and structure. There is a leading acoustic wave which 
propagates outside of the jet toward the farfield and toward the wall, as well as upstream both 
inside and outside the nozzle. As this wave reaches the wall, it reflects back and crosses the shear 
layer, trailing the initial disturbance, as shown in the figure. A large scale structure, related to 
the spatial instability of the jet shear layer is generated along the jet exit. This structure slowly 
propagates downstream. As a result of the interaction between the flexible panels with the pulse, 
the resulting vibration creates an acoustic disturbance which propagates into the radiation domain. 
In the jet flow field, as the process evolves, additional instability is created behind the initial leading 
wave as observed in the proximity of the nozzle exit. The pressure field is clearly asymmetric about 
the jet axis due to the presence of the wall. 

As the time evolves (Figure 3b) additional pressure waves are created from the jet instability 
which propagate into the farfield. Figures 3b-3e show the evolution of the instability wave as it 
propagates downstream. The additional sound generated from this large scale structure can also be 
seen by following the progression of the leading acoustic wave in Figures 3b-3d. While the initial 
waves due to the starter pulse exit the computational domain (Figure 3d), additional instability 
waves continue to be generated within the jet shear layer, together with the associated shedding of 
acoustic disturbances from the nozzle lip. Hence, long after the initial waves generated by the pulse 
starter have propagated far downstream from the region of interest, indeed, after they have exited 
the domain, sound is continually generated which results in a continual excitation of the panels. 
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In addition, Figures 3a-3e show waves propagating between the pipe and wail boundary, upstream 
of the nozzle exit. The behavior described above has also been verified by a detailed animation of 
successive contours generated by the computation. 

In Figure 4 we illustrate the shedding of disturbances from the nozzle lip. The figure shows 
the pressure in a small region around the nozzle lip for three closely spaced values of time. In the 
top illustration, corresponding to the earliest time, the formation of a pressure disturbance at the 
upper lip of the nozzle exit can be seen. The middle illustration, at a slightly later time, shows 
that this disturbance has shed from the upper lip, while a new disturbance is being generated 
at the lower lip. In the bottom illustration, the disturbance has shed from the lower bp, while 
a new disturbance is forming at the upper lip. The figure illustrates the alternating shedding of 
pressure disturbances from both sides of the nozzle. This behavior is consistent with experimental 
observations [12] and is a source of sound forcing the panels to vibrate. A similar phenomenon 
was found in the axisymmetric computations [21]. Note that for the axisymmetric case, ring type 
of disturbances can be shed from the nozzle lip. The figures also illustrate the propagation of 
disturbances up the pipe. 

3b. Responses of the Panels 

We first consider the PSD of the pressure incident on the centers of each panel. In order to clarify 
the mechanisms influencing the spectral content, we compare the results to a computation without 
flow, in which panel excitation occurs from the starter pulse, as well as to a computation with flow 
but without a bounding wall. The data are presented in Figures 5a and 5b for panels A and B 
respectively. 

Figure 5a shows that the PSDs of the incident pressure on panel A for the wall/flow computation 
are comprised of relatively discrete spectral bands. The behavior for the wall/no flow computation 
exhibits similar discrete bands with lower amplitude levels. In contrast, the computation of no 
wall /flow exhibits a continous smooth spectrum. Thus, the discrete spectral bands exhibited for the 
wall/flow and wall/no flow computations appear to be a direct consequence of repeated reflections 
between the wall and the nozzle. Note that the no wall/flow spectrum is similar to that of the 
farfield pressure at 90° (Figure 2). This is a consequence of the spatial location of panel A relative 
to the nozzle exit. 

The PSDs of the pressure incident on panel B with and without flow exhibit much greater 
differences in comparison to those for panel A. In particular, the flow computations exhibit signif- 
icant amplification of low frequencies as compared to the no flow computation. This low frequency 
amplification is due to sound generated by instabilities within the flow. Note that the PSD of the 
no wall/flow case is similar to that of the farfield pressure at 30° (Figure 2), a consequence of the 
fact that the center point of panel B is nearly 39° from the nozzle exit. The wall/flow computa- 
tion exhibits discrete spectral bands for higher frequencies. We interpret this behavior as due to 
repeated reflections between the wall and the jet boundary. The wall/no flow computation exhibits 
a significantly less pronounced band structure, presumably due to the fact that there is now no 
direct mechanism to promote repeated reflections at the wall. 

Since the panels are identical in all respects, these results indicate that the loading on the panels 
depends very crucially on the geometric placing of the panels relative to the nozzle exit. 

The enhanced low frequency loading on panel B , due to instability wave generated sound, is in 
a frequency band for which the panel can be readily excited. As will be seen, this enhanced loading 
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is reflected in enhanced vibration of and radiation from panel B. We note that these results also 
indicate that near the wall rapid variations in x occur in the unsteady pressure near the nozzle exit. 

The effect of the enhanced low frequency excitation of panel B can be seen in Figures 6 and 
7 where we plot the PSD for the panel velocities and transmitted pressure respectively at the 
panel centers. Both panels exhibit relatively narrow, nearly coincident bands of peak spectral 
response. We denote these bands for the velocity spectra by fl v \ , A, 3 and for the transmitted 
pressure spectra by /3 p \, fi p 2, fi p 3. The enhanced low frequency response of panel B is apparent 
from both figures. In view of the analysis of the incident pressure above we attribute this enhanced 
response primarily to instability wave generated sound which provides enhanced excitation of the 
downstream panel B. We note that the broadband incident pressure on panel B is converted to a 
discrete velocity and transmitted pressure spectra. The three primary pressure bands, /3 P 1, f) p 2, /? P 3, 
are clearly associated with the analogous velocity bands, [i v 1 , f) v 2, The large response for panel 

B in the band fi v 2, together with the associated radiation (/i P 2), represent the main effect of the 
instability wave generated sound on the panels and reflects the frequencies within the broadband 
indicated in Figure 5b at which the panel is most responsive. 

The results indicate the role of panels located downstream of the jet exit as a narrow filter, 
converting a broadband forcing into relatively narrow spectral bands. The bands are concentrated 
in the low frequency range, characteristic of instability wave generated sound, which indicates that 
shear layer instabilities can be an important mechanism in panel excitation. 

3c. Acoustic Radiation from the Panels 

The time evolution of the pressure radiated by the motion of the panels is shown at successive 
times in Figures 3a- 3e. Figure 3a shows a nearly circular leading wave radiating away from the 
panels. This is due to the excitation of the panels by the leading acoustic disturbance initiated by 
the pulse starter. In addition, there is an asymmetrical disturbance generated primarily from panel 
B which propagates mainly in the downstream direction (Figure 3b). 

As time evolves, repeated radiation from both panels occurs. The figures clearly demonstrate 
a preferred downstream beaming of this radiation, although there is noticeable upstream acous- 
tic generation as well. This is in contrast to the nearly circular disturbance generated by the 
pulse starter. We note the presence of well defined radiation even at late times, after the initial 
disturbance due to the pulse starter has exited the computational domain. 

PSDs of the radiated pressure at sample farfield points are shown in Figure 8. Results are 
presented along the line y — 36 D for selected values of x. The farfield spectra exhibit relatively 
discrete spectral bands. The three pronounced frequency bands (fi p 1, f) p 2, /i P 3 ) , are analogous to 
those for the velocity and transmitted pressure (Figures 6 and 7). The primary feature of the 
radiation is the presence of the large peak (fi p \) near 250 Hz. This is an effect of the flow and, 
based on the transmitted pressure (for example, Figure 7), is primarily due to radiation from the 
downstream panel (panel B). We note that while the panels can be most readily excited at low 
frequencies, the comparable peak in a computation with no flow is roughly three times smaller than 
the peak indicated in Figure 8. Thus, the low frequency spike is a consequence of low frequency 
sound generation within the flow field. The enhanced radiation corresponding to this peak is 
primarily a manifestation of the excitation of the panels by instability wave generated sound. 
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4. Conclusion 


We have computed the full flow/acoustic/structure coupling for a model of a panel-stringer assembly 
forced by jet noise. The jet is initially excited by a pulse and continues radiating sound long after 
the initial disturbance due to the pulse has exited the computational domain, thus simulating the 
behavior of a real jet flow field. This radiation is due to the triggering of instability waves in the 
jet. The instability waves act as sources of sound, generating directionally beamed low frequency 
sound into the farfield primarily at low to mid angles to the jet axis. Instability wave generated 
sound acts to excite the flexible panels, particularly the panel located downstream of the nozzle 
exit. Thus, the location of the panels with respect to the jet exit is critical in terms of the resulting 
structural vibration and acoustic radiation. 

The panels act as a narrow filter, converting the relatively broadband incident pressure to a 
transmitted pressure exhibiting relatively narrow spectral bands. The narrow spectral bands are 
also present in the far field radiation from the panels. The radiation from the panels is beamed 
primarily downstream and vertically from the panels. The frequencies of the radiated pressure are 
related both to properties of the panels and to the properties of the incident forcing pressure, which 
in turn is dependent on the shear layer instabilities of the jet. Due to a continual excitation of 
instabilities in the jet, radiation persists long after the initial disturbance generated by the pulse 
exits the computational domain. 
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Fig. 1 Computational Domain 


Far Field Boundary 
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Noimalized PSD 


Fig. 2 PSD of farfield p—Poo on a circle of radius r = 20 D 
at 30 °, 60 °, 90 ° in a domain with no bounding walls 
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Fig. 3 Contours of constant p—poo at successive times 

(b )t = 38.4D/COO 
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Fig. 3 Contours of constant p—poo at successive times 

(e)f = 73.3Z>/coo 
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Fig. 4 Evolution of pressure disturbances shed from 
the nozzel at 3 successive times (from top to bottom) 

t = 34.9-D/coo, t = 36.6D/COO, t = 38.4.D/COO. An expanded 
region about the nozzle lip is shown. 
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Fig. 5 PSD of incident p — poo at the center of (a) panel A. 



Frequency (Hz.) 
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Fig. 5 PSD of incident p — poo at the center of (b) panel 
B for the computations of wall/flow, wall/noflow and 
no wall /flow. 



Noimalized PSD 


Fig. 6 PSD of v at panel centers 
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Normalized PSD 


Fig. 7 PSD of transmitted pressure 
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Fig. 8 PSD of radiated p — at 

x = -12 D,x = = 18D, along the line y = 36 D. 
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